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1,  INTRODUCTION 


A  cylinder  of  explosive  with  u  hollow  cuviiy  in  one  end  and  u  detonator  at  the  opposite  end  is  k.nown 
as  a  shaped  charge.  The  hollow  cavity  (which  may  assume  almost  any  geometric  shape  such  as  a 
hemisphere,  cone,  tulip,  trumpet,  or,  in  fact,  any  arcuate  device)  is  usually  lined  with  a  thin  layer  of  metal. 
The  liner  forms  a  Jet  when  the  explosive  charge  is  detonated.  Upon  detonation,  a  spherical  wave 
propagates  outward  from  the  point  of  initiation.  This  high-pressure  shock  wave  moves  at  a  very  high 
velocity,  typically  around  1(1  km/s.  As  the  detonation  wave  engulfs  the  lined  cavity,  the  material  is 
accelerated  under  the  high  detonation  pressure,  collapsing  the  liner.  During  this  process  the  liner  material 
is  driven  to  very  violent  distortions  over  very  short  time  intervals,  at  strain  rates  of  lO**  to  10^  s~^  and 
peak  strains  greater  than  10  are  possible.  The  Jet  temperature  is  about  500  to  900  K  according  to  the 
measurements  of  von  HoUe  and  Trimble  and  there  may  be  temperature  gradients  along  and  through  the 
Jet.  The  collapse  of  the  conical  liner  material  onto  the  centerline  forces  a  portion  of  the  liner  to  flow  in 
the  form  of  a  Jet  where  the  jet  tip  can  travel  in  excess  of  10  km/s.  Because  of  the  presence  of  an  axial 
velocity  gradient,  t.he  Jet  will  stretch  until  it  fractures  into  a  column  of  particles.  This  fracture  of  the  Jet 
into  a  scries  of  particles  is  termed  Jet  breakup  or  particulation. 

The  penetration  of  a  shaped  charge  Jot  into  most  target  materials  increases  to  a  maximum  and  then 
decrease.s  as  the  standoff  distance  (distance  from  the  front  of  the  shaped  charge  to  the  target)  increases. 
This  penetration  peak  occurs  Just  prior  to  the  onset  of  Jet  breakup  due  to  the  dispersion,  spread,  and 
tumbling  of  the  Jet  particles  after  particulation.  As  a  result,  it  would  be  advantageous  to  the  shaped  charge 
designer  to  predict  and  control  the  jet  breakup  time. 

The  nature  of  the  particulation  of  the  shaped  charge  Jet  has  been  of  interest  for  over  50  years.  It  is 
known  that  some  materials,  namely,  certain  metals,  exhibit  extreme  ductility  under  the  Intense  dynamic 
conditions  involved  in  the  shaped  charge  collapse  process.  These  materials  often  do  not  possess  the  sumo 
degree  of  ductility  under  ambient  conditions  and  undergo  dynamic  elongations  of  1,000%  or  more.  The 
problem  is  complicated  by  the  fact  that  the  material  properties  of  the  liner  are  not  well  known  under  the 
intcn.se  dynamic  conditions  that  the  Jet  undergoes  during  its  collapse,  formation,  and  growth.  Complex 
hydmeode  computer  programs  arc  limited  because  accurate  equations  of  state  and  constitutive  equations 
arc  not  available  under  these  conditions.  Also,  the  fracture  mechanism  and  associated  algorithm  is  not 
well  known.  Nevertheless,  shaped  charge  experiments  provide  an  excellent  test  bed  for  the  study  of 
materials  under  intense  loading  conditions. 


I 


2.  EXPERIMENTAL  APPROACH 

The  cxperimcnial  characterization  of  shaped  charge  jets  Is  most  commonly  performed  using  multiple 
flash  x-ray  units.  Each  Hash  x-ray  unit  provides  an  image  of  the  shaped  charge  Jet  at  a  Icnown  time.  In 
a  typical  experiment,  several  x-ray  units  arc  flashed  at  prc-determined  time  intervals.  The  x-ray  films  arc 
then  analyzed  to  determine  the  position,  length,  radius,  mass,  and  velocity  of  each  of  the  Jet  particles.  In 
general,  the  entire  shaped  charge  Jet  cannot  be  captured  in  a  single  experiment  due  to  limitations  in  the 
length  of  x-ray  film  which  con  be  exposed  and  the  need  to  protect  the  film  from  the  explosive  blast.  The 
data  obtained  from  the  x-my  films  is  utilized  to  calculate  the  Jet  breakup  time, 

The  Jet  breakup  time  can  bo  experimentally  determined  in  several  ways  with  varying  results.  'Hie  most 
common  method  of  calculating  an  aggregate  Jet  breakup  time  is  to  divide  the  cumulative  Jet  length  by  the 
change  in  the  Jet  velocity  from  the  tip  particle  to  the  slowest  Jet  particle  characterized  in  the  experiment. 
This  "slowest"  particle  may  or  may  not  be  the  rear  of  the  Jet.  This  method  provides  a  single  value,  termed 
the  "cumulative  breakup  time,"  for  the  entire  Jet.  The  cumulative  breakup  time  Is  based  on  several 
assumptions.  The  jet  is  assumed  to  stretch  from  an  initial  length  of  zero  at  a  constant,  uniform  rate.  In 
addition,  the  Jet  Is  assumed  to  break  simultaneously,  from  tip  to  toil,  at  the  time  it  reaches  its  maximum 
length. 

The  cumulative  breakup  time  provides  an  accurate  measure  of  the  total  Jet  length  available  for 
penetration  if  it  is  observed  that  the  cumulative  Jet  length  is  a  linear  function  of  the  velocity  of  the  Jet 
panicles.  However,  in  many  cases,  the  cumulative  jet  length  is  nonlinear  with  respect  to  the  Jet  velocity. 
In  such  cases,  the  value  of  the  cumulative  breakup  time  will  vary  with  the  velocity  of  the  slowest  Jet 
particle  included  in  the  calculation.  For  example,  Figure  1  is  a  plot  of  the  cumulative  breakup  time  as  a 
function  of  the  velocity  of  the  slowest  Jet  particle  included  in  the  calculation.  This  plot  was  derived  from 
a  single  experiment  involving  the  Ballistic  Research  Laboratory  (BRL)  81 -mm  standard  shaped  charge. 
The  velocity  of  the  tip  particle  was  7,7  km/s,  and  the  velocity  of  the  slowest  particle  characterized  in  the 
experiment  was  2.1  km/s.  The  cumulative  breakup  time  for  this  experiment  was  found  to  be  147.8  .us, 
The  data  plotted  were  determined  us  the  sum  of  the  lengths  of  the  first  n  particles  divided  by  the  velocity 
difference  between  the  first  particle  and  the  n^  particle.  The  axial  velocity  in  Figure  1  reflects  the 
velocity  of  the  n'^  particle.  The  extremely  large  breakup  times  observed  in  the  leading  portion  of  the  Jet 
arc  due  in  part  to  the  deceleration  of  the  lip  particle  which  reduces  the  difference  In  velocity  between  the 
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tip  particle  and  the  particles  immediately  following.  As  evidenced  in  Figure  1 ,  the  measured  value  of  the 
cumulative  breakup  time  is  dependent  on  the  experimental  setup.  If,  for  example,  the  experimental  setup 
were  such  that  only  the  jet  particles  travelling  faster  than  4  km/s  were  characterized,  tlten  the  cumulative 
breakup  time  would  be  reported  as  132  ps. 

A  more  accurate  representation  of  the  experimental  data  may  be  obtained  by  treating  the  cumulative 
Jet  length  as  a  piecewise  linear  function  of  the  Jet  velocity  (Held  1985).  Thus,  the  Jet  breakup  time  may 
be  calculated  for  several  velocity  ranges.  In  fact,  the  breakup  time  can  be  calculated  for  each  particle. 
For  example,  Figure  2  is  a  plot  of  the  individual  breakup  times  calculated  for  the  same  81 -mm  standard 
shaped  charge.  The  individual  breakup  time  is  extremely  sensitive  to  nonuniform  velocity  differences 
between  particles  and  to  measurement  error.  Some  of  the  Jet  particles  separate  from  one  another  at 
breakup  but  travel  at  very  nearly  identical  velocities  after  separation,  thus  allowing  unrealistically  large 
absolute  values  of  breakup  times  to  be  computed.  Breakup  times  which  are  calculated  over  large  or 
predetermined  velocity  intervals  arc  inherently  more  .stable. 

The  breakup  times,  determined  from  analysis  of  the  flash  radiographs,  discussed  so  far  are  based  on 
the  assumption  that  the  Jet  stretches  at  a  constant  rate  until  it  reaches  a  ma.ximum  length  and  particulates. 
The  reference  time  for  this  process  is  the  point  In  time  and  space  from  which  the  Jet  emanates,  known  as 
the  "virtual  origin."  Conversely,  the  breakup  time  can  also  be  calculated  under  the  assumption  that  the 
Jet  particles  travel  at  a  constant  velocity  after  the  Jet  particulates.  Thus,  the  time  at  which  two  Jet  particles 
separate  from  one  another  can  be  derived  from  the  distance  between  the  two  particles  at  a  known  time  and 
the  velocity  of  the  particles.  The  reference  time  in  this  case  is  determined  by  the  experimentalist  (i.e.,  the 
reference  point  may  be  the  activation  of  the  detonator,  the  time  the  detonation  wave  reaches  the  liner  apex, 
the  time  the  jet  tip  reaches  the  base  of  the  liner,  or  other  reference  points).  Figure  3  is  a  plot  of  the 
separation  times  computed  from  the  81 -mm  sfiapcd  charge  experiment.  The  separation  times  plotted  in 
Figure  3  arc  relative  to  the  activation  of  the  detonator.  The  separation  limes  tend  to  increase  with 
decreasing  Jet  particle  velocity,  which  indicates  the  Jet  broke  from  the  tip  to  the  tail,  although  in  a 
nonuniform  manner.  The  separation  time  calcuiation  also  tends  to  give  unrealistic  values  when  the 
difference  in  velocity  between  particles  is  nearly  zero.  However,  unlike  the  individual  breakup  time 
calculation,  the  separation  time  is  measured  as  a  very  small  value  (or  large  negative  value)  as  the  vefxtity 
difference  approaches  zero. 
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Also,  the  distance  a  jet  particle  travcis  before  separation  from  its  adjacent  particies  can  be  measured. 
The  measurement  of  distance  to  particulation  is  very  similar  to  the  separation  time  measurement  and  is 
based  on  the  same  assumptions.  Figure  4  is  a  piot  of  the  separation  distance  vs.  jet  velocity  again  from 
the  same  8 1 -mm  shaped  cliargc  experiment.  The  distance  to  separation  is  shown  to  decrease  with 
decreasing  jet  velocity.  'Ilie  distance  to  separation  is  measured  relative  to  the  liner  base. 

3.  ANALYTICAL  MODELS:  BACKGROUND 

Analytical  models  of  the  jet  breakup  arc  available,  notably  due  to  Walsh,  Held,  Carieone,  Chou, 
Pferfer,  Haugstad,  and  Hirsch,  but  are  primarily  one-dimensional  and  semi-empirical,  More  theoretical 
models  are  available  from  Pack,  Curtis,  Romero,  Frankel,  Welhs,  and  others.  Early  attempts  at 
understanding  the  breakup  of  a  shaped  charge  jet  involved  studies  of  liquid  jet  breakup  since,  visually, 
the  breakup  phenomenon  would  appear  to  have  .some  analogy  to  the  breakup  of  a  liquid  jet  as  studied  by 
Rayleigh  (1894).  More  recent  liquid  jet  breakup  studies  concern  atomization  or  the  breakup  of  liquid 
drops  in  air  where  fluid  vi.scosity  and  surface  tension  ore  the  dominant  mechanisms.  Such  analyses  are 
given,  for  example,  by  Gordon  (1959)  and  Heidmann  and  Grocneweg  (1968).  Tomotika  (1935,  1936) 
considered  the  analysis  of  an  incompressible,  viscous  Newtonian  fluid  surrounded  by  another  viscous  fluid 
under  the  action  of  capillary  (viscous  and  surface  tension)  forces.  A  stability  analysis  was  performed  for 
the  governing  equations  which  accounted  for  stretching  jets  but  neglected  inertia  effects.  Mlkami  et  al. 
(1974)  studied  the  breakup  of  a  viscous  liquid  thread  in  a  viscous  liquid  (pulp  fibers  in  water)  for  the 
paper  pulp  industry.  Again,  a  stability  analysis  was  performed  for  stretching  jets  in  the  absence  of  inertia 
terms.  Goldin  et  al.  (1969)  considered  the  breakup  of  a  viscoelastic  fluid  column,  and  extended  the  above 
analyses  to  a  non-Newtonian  flow.  Frankel  and  Welhs  (1985)  attempted  to  directly  consider  the  breakup 
of  a  shaped  charge  jet  based  on  the  studies  described  above  but  included  inenial  effects.  The  shaped 
charge  jet  was  assumed  to  be  liquid  with  a  surface  tension  of  3  x  10^  N/m.  A  perturbation  solution  was 
performed,  but  it  was  ob.scrvcd  that  the  jet  must  have  a  very  high  yield  strength  (in  excess  of  2.4  GPa  for 
copper)  for  the  surface  tension  in  the  analysis  to  be  meaningful.  He  noted  that  such  high  values  of  the 
dynamic  yield  .strength  of  copper  were  quoted  by  Van  Thiel  and  Levatin  (1980).  Curtis  (1987)  noted  that 
the  surface  tension  required  to  produce  significant  instability  growth  was  of  the  order  of  the  jet  yield 
strength  divided  by  the  jet  radius  which  is  several  orders  of  magnitude  (namely,  4  x  10^)  greater  than  the 
known  surface  tension  of  copper  in  air  at  ambient  conditions.  Therefore,  the  tensile  stress  in  the  jet  was 
considered  to  be  the  mechanism  in  the  jet  which  produced  the  instability  growth.  Curtis  used  tlie 
equations  of  motion,  the  irrotatlonality  condition,  and  the  von  Mises  flow  criterion,  and  he  expressed  the 
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stress  tensor  as  a  superposition  of  a  pressure  term  and  a  uniaxial  stress  term.  A  perturbation  and  stability 
analysis  was  performed  for  an  assumed  (zeroth  order)  velocity  field.  It  was  noted  that  the  Jet  breakup  is 
a  function  of  the  jet  strength  to  density  ratio,  but  the  compicxity  of  the  equations  precluded  a  detailed 
expression  for  the  jet  breakup  time.  Curtis  concluded  that  it  would  be  more  appropriate  to  use  the  Levy- 
Miscs  flow  criterion  (the  full  classical  plastic  theory)  instead  of  the  von  Mises  criterion.  Pack  (1987) 
analyzed  the  plastic  How  problem  by  using  the  equations  of  motion  and  the  Levy-Mises  stress-strain 
relationships  plus  the  von  Mises  criterion  for  plastic  flow.  A  perturbation  analysis  using  a  scries  solution 
was  conducted  along  with  a  stability  analysis,  again  for  an  assumed  velocity  field.  Pack  then  postulated 
two  causes  of  breakup:  1)  necking  of  the  Jet,  and  2)  the  creation  of  voids  in  the  Jet  when  all  principal 
stresses  become  tensile.  His  predicted  results  for  the  two  causes  were  very  close  and  he  was  not  able  to 
distinguish  between  the  two  mechanisms.  Romero  (19S9),  along  the  same  lines,  performed  a  stability 
analysis  of  a  rapidly  stretching  Jet  and  concluded  that  the  Jet  stability  depends  on  a  parameter  he  called 
"F,"  which  represents  the  ratio  of  the  inertial  forces  to  the  plastic  forces.  In  fact,  P  is  time  dependent  and 

r2(t)  -  J^pp2(t)a2(t), 
o 

where  a  Is  the  Jet  radius,  p  is  the  strain  rate,  p  is  the  Jet  density,  and  a  is  the  yield  stress,  The  Jet  is  stable 
until  r  is  less  than  1.  Romero's  analysis  was  much  like  that  of  Frankel  and  Welh's,  but  used  the  Levy-. 
Mises  equations  for  a  perfectly  plastic  material.  As  in  the  earlier  analysis,  a  perturbation  solution  was 
obtained. 

Other  studies  of  the  Jet  breakup  time  were  conducted  by  Shelton  and  Arbucklc  (1979)  who  considered 
the  propagation  of  relief  waves  following  a  break  in  the  Jet.  The  speed  of  relief  wave  propagation  from 
the  break  was  calculated  by  two  different  models.  Walsh  (1984)  provided  insight  into  the  mechanisms 
(I.C.,  surface  disturbances)  that  may  lead  to  Jet  instability.  These  disturbances  include  explosive 
homogeneity,  liner  dimensional  tolerance,  velocity  gradient  perturbations,  and  liner  dynamic  strength 
variations.  Walsh  concluded  that  breakup  depends  on  the  perturbation  structure  and  a  single  dimensionless 
flow  parameter, 

A  - 

0  ■ - f-, 

AV'R 
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where  o  is  the  yield  strength  of  the  jet,  p  is  its  density,  R  is  the  jet  radius,  and  is  the  initial  jet 
stretch  rate.  The  subscript  o  represents  initial  jet  condidons.  Finally,  it  was  shown  that  the  length  traveled 
until  breakup  occurred,  L),,  is  given  by 

where  Lg  is  the  initial  jet  length,  (|)q  Is  the  initial  value  of  <(>,  and  'Fq  is  the  initial  surface  roughness. 

Analytical,  one-dimensional  models  arc  also  available.  Miller  (1982)  showed  that  a  one-dimensional 
theory  applied  to  a  perfectly  plastic  stretching  metal  rod,  or  applied  to  a  stretching  metal  rod  governed  by 
the  Steinberg-Ouinan-Cochran  constitudve  model,  predicts  the  generation  of  a  progression  of  new  necks 
from  an  existing  neck.  Two-dimensional  finite-difference  calculations  predict  the  same  behavior. 

Chou,  Sidhi,  and  Mortimer  (1963),  in  perhaps  one  of  the  earliest  analytical  studies  of  shaped  charge 
jet  breakup,  deduced  three  possible  mechanisms,  or  model  approaches,  for  breakup.  These  mechanisms 
were  viscoplastic  effects  (which  as  Chou  et  al.  noted,  would  preclude  Hopklnson  scaling  of  the  jet 
breakup— if  this  is  a  limitation);  a  Mott-type  statistical  approach  which  assumes  that  the  tensile  strength 
is  not  constant  but  exhibits  a  scatter  characteristic  (l.e.,  there  exists  a  scatter  in  the  value  of  the  reduction 
in  area  (or  value  of  strain)  at  which  fracture  occurs  in  a  tensile  test);  and  finally,  temperature  or  molten 
jet  effects,  Chou  et  al.  noted  that  any  or  all  of  these  mechanisms  could  influence  jet  breakup.  In  later 
publications,  Chou  and  Carlcone  (1977b)  developed  a  one-dimensional  Lagrangian  theory  for  stretching 
plastic  jets  with  a  constitutive  equation  of  the  form 


o  -  Og  +  Ce  +  Re^  + 


where  C,  R,  etc.,  are  positive,  empirical  coefficients.  The  stretching  jet  problem  assumed  small  surface 
disturbances  and  small  strains.  A  stability  analysis  was  performed  using  a  linearized  form  of  the  one¬ 
dimensional  model  and  a  stability  criterion  was  established.  The  theory  showed  that  the  disturbance 
growth  rate  in  the  jet  Increases  as  the  wavelength  of  the  disturbance  decreases  and  as  the  stretching  rate 
decreases.  Also,  for  a  perfectly  plastic  jet  material,  where  o  -  Og,  the  disturbance  growth  rate  increases 
with  increasing  o/p.  All  trends  predicted  by  the  one-dimensional  theory  were  in  agreement  with  two- 
dimensional  hydrocode  calculations.  Also,  the  hydrocode  calculations  revealed  the  existence  of  a  critical 
wavelength  (l.c.,  disturbances  having  this  wavelength  will  grow  faster  than  all  others)  and  showed  that 
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elasticity  effects  and  compressibility  effects  were  nearly  negligible.  Additional  details  and  analyses  are 
given  in  Chou,  Carleone  and  Karpp  (1974)  and  Chou  and  Carleone  (1976).  Chou  and  Caileone  (1977a, 
1977b)  provide  details  regatding  the  Jet  growth  function  and  other  factors  affecting  stability.  Qiou  and 
his  coauthors  approached  the  calculation  of  the  breakup  time  directly  utilizing  many  Hash  radiographs  of 
jets  from  shaped  charge  liners.  The  final  results  of  the  (Thou -Carleone  model  for  a  copper  jet  reveal 


Tb  -  3.75  -  0.125  to  + 

Eo 


where  Tb  «  dimensionless  breakup  time  «  Cp  Cp  -  ^o/p  ;  Jo  ■  dimensionless  strain  rate  ■ 
;  and  t^  »  initial  strain  rate  »  AV/AX;  and  ^  is  the  breakup  time,  a  is  the  yield  strength,  AV 

is  the  change  in  velocity,  and  AX  is  the  change  in  length.  The  subscript  o  designates  initial  values  when 
the  liner  clement  first  arrives  at  the  axis  of  symmetry.  The  numerical  coefficients  result  from  a  curve  fit 
of  the  dynamic  ductility  factor  to  experimental  data.  The  details  of  this  model  and  comparison  with 
experimental  data  are  summarized  in  Walters  and  Zukas  (1989). 

Recently,  Chou  et  al.  (1992)  expressed  the  jet  breakup  time  as 


f  2  ^ 

1/3 

(  2  ^ 

1/3 

f  ■'1 

■  k 

l-olo 

■  k 

r^t 

kJ 

1/3 


where  the  radius  r  and  time  t  of  a  jet  segment  may  be  taken  at  any  time  before  jet  breakup.  The  factor 
k  Is  taken  to  be  5.0  based  on  a  fit  of  the  equation  with  experimental  data.  The  terni  Cp  is  defined  as 
before,  but  o  is  the  yield  stress  (Oy)  for  an  isotropic  jet  and  Oy  (0.5  -t-  1/m^)  for  an  anisotropic  jet,  where 
m  is  the  ratio  of  transverse  strength  to  axial  strength.  Also,  Rottenkoller(1989)  obtained  a  similar  formula 
for  the  breakup  time  of  a  shaped  charge  jet. 


\/> 
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where  dg  is  the  initial  jet  diameter,  and  B  is  equal  to  ^d/Cp  and  assumed  to  be  constant. 


Eitan  Hirsch  (1979)  expressed  the  breakup  time  as 


X 


djo/Vp,. 


where  djo  is  the  Initial  Jet  diameter  when  elongation  staits,  R  is  the  radius  of  the  liner  element,  t  is  the 
thickness  of  the  element,  p  is  the  collapse  angle,  x  is  the  breakup  time  measured  from  the  arrival  of  the 
explosive  wave  front  at  the  liner  element  from  where  the  jet  originates,  and  Vp|  is  a  velocity  specified  by 
the  liner  material,  called  the  "plastic  velocity."  Hirsch  (1979)  expressed  the  plastic  velocity  as  the  velocity 
difference  between  the  Jet  particles  or  Vp|  •  (Vj  -  V^)  /  n,  where  n  Is  the  number  of  jet  particles,  Vj  is  the 
jet  tip  velocity,  and  V,  is  the  velocity  of  the  rear  of  the  jet.  Hirsch  also  gave  expressions  for  the  breakup 
time  of  an  expanding  cylinder  (pipe)  and  a  linear  shaped  charge.  Hirsch  (1989)  attempted  to  express  the 
plastic  velocity  us 

Vpi  ■  v^o/p  . 

where  o  is  the  liner  metal  dynamic  yield  stress  and  p  is  the  density  of  the  liner  material,  from  microscopic 
metallurgical  conditions.  Hirsch  employed  the  Mott  fragmentation  model  as  mentioned  by  Chou  et  al. 
(1963)  and  provided  a  description  of  the  breakup  time  model  as  arising  from  shear  bands  or 

Vpi  -  /do^P. 

where  do^  represents  the  difference  between  the  isothennal  and  adiabatic  stress  vs,  strain  characteristics 
of  the  metal  at  the  point  where  adiabatic  stress  becomes  a  maximum.  Hirsch  quotes  strains  of  1  to  2  ut 
a  strain  rate  of  U)^  s~’  for  OPHC  copper.  These  values  of  strain  are  considerably  lower  than  those  quoted 
by  Chou  and  Carieone  (1976).  Hirsch  (1981a)  hirther  qualified  the  plastic  velocity  by  suggesting  a 
breakup  mechanism,  where  holes  caused  by  a  pile  up  of  vacancies  are  formed  at  the  metal  .surface  and 
gradually  increase  until  breaking  is  caused  by  the  formation  of  voids  in  the  jet.  Hirsch  also  predicted  the 
existence  of  a  strain  rate  threshold  below  which  other  mechanisms  dominate  the  breakup  process.  Hirsch 
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(1981b)  states  that  even  perfectly  symmetrical  and  homogeneous  shaped  charge  conflgurations  have 
transverse  velocity  components  in  the  Jet.  This  means  that  the  breakup  process  starts  during  the  liner 
collapse  and  the  transverse  velocity  influences  the  Jet  breakup.  Hirsch  (1990)  relates  the  plastic  velocity 
to  the  proccs.ses  which  affect  the  liner  metallurgical  state  during  the  initial  stoges  of  jet  formation.  Hirsch 
shows  how  both  the  deformation  energy  heating  the  sliding  shear  bands  during  the  localization  process 
and  the  rate  of  the  instability  growing  in  the  plastic  flow  during  this  process,  combine  to  determine  the 
plastic  velocity  parameter.  This  velocity  is  shown  to  be  related  to  both  the  velocity  due  to  the  plastic 
deformation  and  the  component  of  the  maximum  slide  velocity  allowable  to  form  shear  bands  in  the 
elongation  direction.  Hirsch  thus  attempts  to  Include  the  influence  of  Ute  metallurgical  structure  of  the 
liner  on  the  breakup  time.  In  this  study,  small  strains  (less  than  1)  are  calculated. 

Haugstad  (1983)  and  Haugstod  and  Dullum  (1983)  formulated  two  models  for  the  jet  breakup  time 
by  considering  vlscoplastlc  effects  or  <j  «  +  pft,  where  p  has  the  form  of  a  viscosity,  a  is  the  yield 

stress,  Oq  Is  the  quasistatlc  yield  stress,  and  i  is  the  strain  rate.  They  obtained 

,  .  1 
t(j  __  “■  —  < 

/Oq/p 

in  the  limit  as  p  -»  0,  where  is  the  initial  jet  diameter  and  a  is  a  constant.  Note  that  if 


a  modified  Hirsch  model  results.  The  second  (viscoplastic)  model  gave 

tb  -  b  p/Oo  -  1/to. 

where  b  is  a  constant. 

pfeffer  (1980)  obtained  u  jet  breakup  time  formula  which  Indicates  the  breakup  time  is  inversely 
proportional  to  the  initial  strain  rate,  weakly  dependent  on  the  Initial  jet  radiu.s,  but  independent  of  the 
liner  material  yield  strength.  Pfeffer  gives 


13 


tb  -  1.4/feo  ♦  48.5  r^/C^, 


where  is  the  initial  strain  rate,  r^  is  the  initial  jet  radius,  and  Cq  is  the  shock  velocity  in  the  jet.  PfeiTer 

assumed  a  formuia  for  the  shape  of  the  broken  jet  segment  and  Ids  results  are  based  on  a  curve  fit  of  two- 
dimensional  computer  simulations. 


Held  (1989)  and  Mayseless  ct  ol.  (1989)  (os  well  os  the  Soviet  researchers)  advocate  a  calculation  of 
the  breakup  distance  instead  of  the  breakup  time.  Walters  and  Summers  (1992a)  derived  the  fonniila  of 
Mayseless  et  al.  (1989)  from  the  velocity  difference  expressions  shown  in  Walters  and  Summers  (1992b) 
and  assuming  simple  tension.  Held,  Mostert  and  Koenig.  Oolaskl  and  Duffy,  and  Hirsch,  among  others, 
comment  on  the  influence  of  liner  metallurgy  on  the  jet  breakup  time.  Oolaski  and  Duffy  (1987)  did  not 
provide  a  breakup  time  formula,  but  showed  a  direct  correlation  between  liner  grain  size  and  jet  breakup 
time.  Mostert  and  Koenig  (1987)  claimed  that  the  jet  from  a  shaped  charge  elongates  to  a  strain  well  in 
excess  of  10  before  it  particulates.  They  also  noted  that  the  micromechanical  properties  of  the  liner,  as 
well  as  its  purity,  have  an  Influence  on  the  ductility  of  the  jet,  but  these  factors  are  not  included  in  the 
existing  breakup  models.  Mostert  and  Koenig  expressed  the  breakup  time,  x,  as 


(t  *  -  «o/Vo)^^ 


ln(No/N,). 


where  the  jet  has  an  initial  diameter  d^,  and  an  initial  length  represents  the  velocity  difference 

ucms.s  the  initial  length  and  a  is  a  proportionality  constant.  is  the  Initial  moving-dislocation  density 
immediately  after  jet  formation  and  is  influenced  directly  by  the  Initial  detonation  shock  pressure  on  the 
liner,  the  pressure  experienced  in  moving  through  the  collapse  zone,  and  the  effective  pulse  duration. 
Al.so,  Nq  varies  inversely  with  the  liner  grain  size,  N,  is  the  postulated  critical  value  of  the  moving- 
dislocation  density  where  further  elongation  is  inhibited.  The  value  for  the  minimum  density  for  plastic 
elongation  could  conceivably  be  linked  to  the  value  of  the  dynamic  equilibrium  density  just  prior  to 
breakup. 

However,  all  of  these  models  arc  cither  incomplete  or  involve  unknown  material  properties  which  are 
usually  a.ssumcd  to  be  constant,  such  us  flow  stress,  dynamic  viscosity  or  the  like.  Also,  some  models 
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involve  unknown  niting  patameters  or  vuriubics.  In  addition,  models,  or  observations,  that  seem 
acceptable  for  one  liner  geometry  and  material  fail  when  applied  to  other  geometries  and/or  liner  materials. 
Thus,  the  various  analytical  models  do  not  agree  with  each  other  or  with  the  experimental  data,  and,  in 
fact,  the  method  of  calculating  Jet  breakup  time  distribution  from  the  experimentally  obtained  na.sh 
radiographs  varies  from  institution  to  institution.  Nevertheless,  significant  insight  has  been  gleaned  from 
the  analytical  breakup  time  models,  namely,  the  existence  of  a  critical  wave  length,  the  effect  of  jet 
strength  (o)  and  density  (p),  and  the  dependence  on  strain  rate. 

It  is  known  experimentally  that  the  jet  microstructurc,  notably  the  grain  si/,c  (and  probably  the  grain 
size  distribution,  the  grain  orientation,  the  metal  chemistry  or  purity,  and  the  material  texture)  strongly 
affects  the  jet  breakup  time.  In  fact,  Chokshi  and  Meyers  (1990)  point  out  that  high  strain  rate 
deformation  leads  to  on  increase  in  temperature  which,  in  conjunction  with  the  large  strains  involved,  leads 
to  a  very  fine  grain  size  microstructurc  due  to  dynamic  rccrystallization.  Subsequently,  the  fine  grain  size 
leads  to  supcrpla,sticity  at  high  strain  rates,  which  in  turn  leads  to  large  tensile  strains  to  failure.  However, 
these  factors  arc  not  included  in  the  analytical  or  theoretical  breakup  time  models  and  including  them  is 
not  straightforward,  in  fact,  the  underlying  mcchanl.sms  behind  Jet  breakup  arc  not  understood.  If  u 
known  perturbation  is  applied  to  a  stretching  Jet  in  a  numerical  experiment,  the  Jet  will  neck  at  a  critical 
wavelength  and  break,  The  exact  nature  and  origin  of  this  critical  perturbation(s)  is  unknown,  but  may 
be  due  to  irregularities  In  the  liner  material  yield  strength  or  other  material  or  micromechanical  properties, 
nonuniformitics  in  the  Initial  Jet  velocity  gradient.  Jet  surface  roughness,  or  due  to  Inherent  perturbations 
In  the  fabiication  of  the  shaped  charge  (c.g.,  materia!  anisotropies,  liner  wall  thickness  variations,  the 
quality  of  the  liner  inner  surface  and  exterior  surface,  inhomogencltics  in  the  explosive  fill,  etc.).  Also, 
the  exact  nature  of  the  fracture  process  (at  a  critical  value  of  strain,  prc,ssurc,  stress,  plastic  work,  or 
internal  energy)  is  unknown.  Furthermore,  models  that  account  explicitly  for  the  nucleation  and  growth 
of  voids,  cracks,  and  shear  bands  have  not  yet  matured  to  the  extent  that  they  con  be  readily  incotporated 
into  the  hydrocode  models. 

Factors  which  are  known  to  affect  Jet  breakup  time  arc  given  by  Held  (1989)  and  Walters  and  Zukas 
(1989).  In  general,  the  Jet  breakup  time  cun  be  increased  by  decreasing  the  Jet  stretching  rate,  increasing 
the  Jet  radius.  Increasing  the  Jet  density,  decreasing  the  Jet  strength,  and  increasing  the  ductility  of  the  Jet 
under  the  dynamic  conditions  described  above.  Thus,  the  liner  design,  liner  geometry,  liner  material,  and 
method  of  fabrication  of  the  shaped  charge  liner  arc  all  pertinent  factors  to  be  considered  in  assessing  the 
Jet  breakup  time  distiibuiion. 
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4.  THE  ANALYTICAL  MODEL 

In  this  section,  un  analytical  model  is  developed  for  the  breakup  time  of  the  jet  from  a  shaped  charge 
iiner.  The  model  is  based  on  three  presumptions.  First,  a  kinematic  expression  for  the  breakup  time: 
second  on  expression  related  to  plastic  stability;  and  finally,  a  mateiiul-bascd  constitutive  equation  relating 
the  stress,  strain,  strain  rate,  and  temperature.  In  other  words,  the  Jet  from  a  shaped  charge  liner  will 
particulate  when  it  becomes  plastically  unstable  (Chou  and  Carleone  1977b),  and  the  breakup  time  will 
depend  on  the  stress,  strain,  strain  rate,  and  temperature  at  failure  (particulation). 

Figure  5  illustrates  the  kinematic  expression  for  the  jet  breakup  time.  An  initial  length  of  jet, 
eventually  stretches  to  length  ( where  it  begins  to  neck  at  the  breakup  time,  x,  Then 


Vr),  (1) 

where  Vj^  is  the  tip  velocity  of  the  jet  and  Vr  is  the  rear  or  tail  velocity  of  the  jet.  Dividing  Equation  1 
by  yields 


t-L 


or. 


X  ■ 


(2) 


The  plastic  stability  criterion  requires  that 


for  stability,  othcrwi.se  the  jet  necks  and  eventually  breaks  (Wallers  and  Summers  1992u,  1992b).  In  order 
to  implement  the  model,  Equation  3  is  used  with  the  equality  sign,  or 
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Figuie5.  The  tinemalic  ctpression  for  the  ii 


~-de.  (4) 

o 

This  is  the  same  equation  used  to  calculate  the  maximum  load  point.  Equation  4  is  integrated  between 
the  limits  e  s  o  to  e  =  ep  and  o  s  to  o  a  Op,  where  Og  is  the  value  of  o  at  e  «  0  and  the  subscript  F 
denotes  values  at  failure.  Thus, 


Op 

_  ■  exp  Sp.  (5) 


Finally,  a  constitutive  equation  which  relates  stress,  strain,  strain  rate,  and  temperature  and  is  valid 
for  high  strain,  strain  rate,  and  temperature  values  is  necessary.  Such  an  equation  may  not  exist  for  the 
extreme  strains,  strain  rates,  and  temperatures  the  Jet  from  a  shaped  charge  liner  is  subjected  to. 
Nonetheless,  two  constitutive  equations  were  considered,  Johnson-Cook  (Johnson  19B3;  Johnson  and  Cook 
1983,  1985)  and  Zerilli-Armstrong  (1987),  due  to  their  popularity  in  many  hydrocodes  and  the  fact  that 
tile  material  coeffleients  are  available  for  most  materials  of  interest  in  shaped  charge  liner  design. 

The  Zcrilli-Annstrong  (1987)  constitutive  equation  for  face-centercd-cubic  materials  is 

o  -  Cg  +  C2e‘'2gxp(-C3T  +  C4Tlnfc).  (6) 

The  stress  o  is  in  MPa,  the  strain  rate  i  is  in  s"',  and  the  temperature  T  is  in  K.  For  an  OFHC  copper 
jet,  Zerilli-Armstrong  gives; 

k  =  microstructural  stress  intensity  =  5.0  MPa  (mm)'^ 

K  =  grain  size  =  0.075  mm 
Cg  =  46.5  MPa 
Cj  =  890.0  MPa 
C3  =  0.0028  K"' 

C4  =  0.000115  K-‘, 
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i'hc  Johnson-Cook  (1983,  1985)  constitutive  equation,  which  has  a  different  functional  form  as 
compared  to  the  ZcriUi-Amistrong  equation  is 

o  -  (A  +  Be^Xl  +  Clnfe)(l  -  T*"’),  (7) 

and 

T  •  -  (T  -  -  T«om)  - 

0  S  T*  i  1, 

where  for  OFHC  copper, 

-  1.356  K.  -  293  K,  C  -  0.025. 

m  ■  1.09,  n  «  0.31,  A  -  90.0  MPa,  and  B  ■  292.0  MPa. 

Note  that  in  the  system  of  units  used  in  the  above  equations,  the  density,  p,  of  copper  is  8,960  kg/m^ 
In  addition,  the  modified  Johnson<Cook  model  was  examined.  The  modified  Johnson<Cook  equation  Is 


o  -  (A  +  Be'')e'(l  -  T**"),  (8) 

where,  for  OFHC  copper,  all  constants  have  the  same  value  as  in  the  original  Johnson-Cook  model,  but 
the  functional  form  of  the  strain  rate  dependency  has  been  changed.  Any  constitutive  equation  may  be 
u.scd  in  this  analysis,  assuming  of  course  that  the  parameters  used  in  the  equation  are  available  and  do  not 
introduce  additional  unknowns  (such  as  the  pressure  for  example). 

The  true  strain  at  failure,  ep,  is  detennined  using  Equation  5  and  a  constitutive  equation  (e.g.. 
Equations  6-8).  The  value  of  is  calculated  as: 

-  C„  +  kX-*^ 
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from  Equation  6  for  the  ZeriilUArmstrong  constitutive  equation  for  face^ccntcred-cubic  materials.  For  the 
Johnson-Cook  constitutive  model,  Equatiuh  7. 

Oo  -  A(1  +  Clnfe)(l  - 

In  addition,  .he  true  strain  rate,  required  in  the  constitutive  equations,  is  detemnined  as  a  hinction  of  the 
true  strain  using 


e  «  exp(-e) 

as  derived  in  Walters  and  Summers  (1992b)  under  the  assumption  of  a  constant  Jet  stretching  rate. 
Finally,  the  jet  temocraiuro  is  assumed  to  be  established  at  failure  and  is  constant.  Thus,  the  tlnal  true 
strain  ut  failure  and  the  true  stress  at  failure  may  be  determined  given  an  Initial  strain  rate  and  an  assumed 
jet  temperature. 

The  initial  strain  rate  is  given  as 


to 


AU 

T’ 


where 


AU  -  V.„  - 


and 

Vjo  =  the  jet  tip  velocity, 

■  the  velocity  of  the  rear  of  the  Jet. 


The  velocily  dirfcrcncc,  AU,  can  be  calculated  using  the  methods  given  in  Walters  and  Zukas  (1989),  for 
example,  and  is  considered  known.  The  initial  Jet  length,  is  taken  to  be  the  initial  length  of  the  Jet 
available  prior  to  stretching.  For  a  conical  liner,  is  determined  based  on  the  slant  length  of  the  cone. 
This  is  essentially  the  Jet  length  considered  in  the  steady  state  Jet  formation  theory  of  Birkhoff, 
MacDougall,  Pugh,  and  Taylor  described  in  Walters  and  Zukas  (1989).  For  nonconical  liners,  is  based 
on  one-half  the  perimeter  of  the  Inside  surface  of  the  liner  obtained  by  taking  a  cross  section  along  the 
axis  of  the  liner.  In  general,  the  initial  jet  length  is  taken  as; 


where  is  the  slant  height  or  one-half  the  liner  perimeter  and  is  the  measured  tail  velocity  from  the 
(lash  x-ray  of  the  Jet  (which  depends  on  the  length  of  Jet  captured  on  the  film).  Thus, 


^0 


The  ratio  (Vj^,  -  V„„)/Vjo  is  introduced  to  account  for  the  fact  that  the  flash  radiograph  may  not  provide 
information  on  the  entire  length  of  the  shaped  charge  Jet.  However,  this  ratio  does  not  account  for  the 
portion  of  the  apex  region  of  a  conical  liner  which  only  contributes  to  the  mass  of  the  Jet  tip  particle  or 
the  portion  of  the  base  region  of  the  liner  which  does  not  contribute  to  the  jet. 

The  John.son-Cook  constitutive  equation  yielded  very  low  values  for  the  strain  at  failure  (as  compared 
to  the  Zcrilli-Armstrong  constitutive  equation).  Also,  because  of  the  functional  fomi  of  the  Johnson-Cook 
equation,  the  strain  at  failure  is  independent  of  the  temperature  and  strain  rate.  The  same  statements  apply 
to  the  modified  Johnson-Cook  constitutive  equation  (Equation  8).  Thus,  the  Johnsoii-Cook  based  results 
will  not  be  presented.  The  tow  strains  obtained  by  the  Johnson-Cook  model  were  also  observed  by  and 
reported  in  Walters  and  Summers  (1992b). 

Figure  6  plots  the  final  true  strain,  ep,  as  a  function  of  the  initial  .strain  rate,  ,  with  temperature  as 
a  parameter.  The  experimental  data  points  are  plotted  for  .several  different  conical  liner  designs  as  well 
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Rnal  True  Strain 


T-1000K 


Figure  6.  The  Hnal  true  strain  vg.  the  Initial  strain  rate  with  temperature  as  a  parameter  for 
several  copper  shaped  charge  liner  designs. 
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as  a  tnimpet  and  a  hemispherical  liner  conflguration.  In  all  cases,  the  liner  material  was  copper.  The 
experimental  points  were  deteimined  based  on  the  calculated  initial  length,  and  the  total  Jet  length, 
L,  measured  directly  From  the  flash  radiographs  of  the  jet.  Note  that  for  the  conical  and  trumpet  liners 
the  average  finai  true  strain  of  the  jet  is  nearly  constant  or 


ep  ■  in 


L 

X 


2.3. 


For  an  average  final  true  strain  of  2.3, 


exp  e  " 


Ao 

T 


f 


J 


for  cylindrical  jet  particles,  or  xlt^  ■  0.32. 

For  the  hemispherical  liner,  the  average  final  strain  is  higher,  about  2.S.  This  difference  In  final  true 
strain  values  may  well  be  due  to  the  different  jet  collapse  and  formation  mechanisms  between 
hemispherical  and  conical  (or  conleal-like)  liners  or  it  may  relate  to  the  method  used  to  calculate  (as 
one>half  of  the  perimeter  of  the  inside  surface  of  the  liner)  for  hemispherical  liners.  Note  that  this  plot 
would  seem  to  Indicate  that  the  jet  temperature  would  be  around  450  K  to  500  K  for  the  conical  liners 
and  somewhat  lower  For  hemispherical  linen  based  on  the  constitutive  equation  used  for  copper.  Note 
also  that  for  any  given  temperature,  the  Final  true  strain  is  nearly  constant  with  strain  rate. 

Table  1  presents  the  experimentally  determined  jet  tip  velocity,  jet  tall  velocity,  total  jet  length,  and 
cumulative  jet  breakup  time  for  several  conical  liner  designs,  a  trumpet  liner  design,  and  a  hemispherical 
liner  design.  The  calculated  total  length  and  cumulative  breakup  time  are  also  shown  for  comparison. 
The  results  presented  in  Table  1  show  excellent  agreement  with  the  experimental  data  at  a  temperature 
of  450  K.  The  experimental  breakup  time  is  derived  from  the  virtual  origin  approximation,  which  gives 
( a  <0  at  the  virtual  origin  time  tg.  Thus,  assuming  tg  ■  0, 
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Then,  if  vk  0,  it  follows  that 


% 


e,  +  1 


(9) 


The  calculated  values  given  in  Table  I  were  computed  using  Equation  9. 


Table  1.  A  Comparison  Between  the  Experimental  and  Analytical  Values  of  Total  Jet  Length 
and  Cumulative  Jet  Breakup  Time  fbr  Several  Shaped  Charge  Liner  Designs 


63.3-mm 

76.2-mm 

81-mm 

101.6-mm 

140-mm 

163-mm 

Cone 

HemI 

Cone 

Cone 

Cone 

Trumpet 

Vtip  (km/s) 

9.2 

4.2 

7.7 

8.8 

8.5 

■Ea 

Vtail 

4.2 

2.4 

2.1 

4.3 

2.6 

Lp  (mm) 

300.1 

279.9 

820.3 

772.1 

1,330.7 

Tg  (ps) 

99.2 

160.9 

147.8 

169.3 

231.3 

Predictions  T  «  430  K  | 

Lj  (mm) 

530.4 

260.1 

WM 

■n 

1,341.7 

Tg  (ps) 

103.2 

149.3 

■isi 

231.7 

Next,  the  cumulative  Jet  breakup  time  Is  plotted  as  a  Ainctlon  of  charge  diameter  for  various 
temperatures  in  Figure  7.  This  plot  was  constructed  by  holding  the  jet  tip  and  tail  velocities  constant  and 
varying  the  liner  diameter.  In  this  case,  the  experimental  results  from  the  81 -mm  conical  liner  were  used. 
The  measured  average  breakup  time  was  147.8  ps  for  this  round,  as  plotted  on  Figure  7.  The  dashed  line 
represents  the  breakup  time  which  would  be  attained  if  homologous  scaling  is  assumed.  In  fact, 
homologous  scaling  of  the  jet  breakup  time  has  been  experimentally  observed  to  be  a  good  assumption, 
usually  accurate  to  within  ±20%.  Figure  7  shows  that  scaling  appears  to  exist  at  temperature  of  400  K 
to  300  K  more  so  than  at  higher  temperatures  (i.e.,  the  slope  of  the  scaling  [dashed]  line  matches  the  slope 
of  the  400-K  line). 


Figure  8  shows  the  predicted  variation  of  the  cumulative  jet  breakup  time  with  temperature  for  a 
81 -mm  diameter  shaped  charge  with  a  coppei  conical  liner.  Note  the  ductility  of  the  jet  is  expected  to 
decrease  with  increasing  Jet  temperature.  This  i.s  due  to  the  constitutive  equation  and  the  observation  that, 
for  facc-centcrcd-cubic  materials,  the  strain  at  maximum  load  decreases  with  increasing  temperature  (Zeriili 
and  Armstrong  1987). 

The  calculation  of  and  dg,  based  on  the  slant  height  of  the  liner  and  the  tip  to  tail  velocity 

difference,  allows  a  quick  calculation  of  the  total  jet  length  and  cumulative  breakup  time.  However,  this 

method  only  provides  a  prediction  of  the  average  strain  attained  by  the  entire  jet  and  a  single  cumulative 

breakup  time.  Alternately,  this  model  may  be  employed  in  conjunction  with  a  Jet  formation  code.  In  the 

jet  formation  code,  the  value  of  is  defined  to  be  the  length  of  a  segment  of  the  Jet  at  the  time  it  reaches 

» 

the  liner  axis  of  symmetry.  The  velocity  difference,  used  to  calculate  initial  strain  rate,  is  defined  as  the 
velocity  of  the  proceeding  jet  segment  less  the  velocity  of  the  succeeding  element.  Ttie  use  of  the  breakup 
time  model  in  conjunction  with  a  Jet  formation  model  also  allows  for  the  prediction  of  individual  breakup 
times,  actual  separation  times,  and  the  distance  a  .segment  of  jet  travels  before  particulation. 

Figures  9-12  repeat  the  experimental  data  shown  in  Figures  l-^  for  the  81 -mm  copper  conical  shaped 
charge  liner.  The  analytical  predictions,  using  the  Jet  formation  code,  are  also  shown  on  these  plots  for 
a  temperature  of  450  K,  Figure  9  is  a  plot  of  the  ctmulatlve  breakup  time  vs.  the  Jet  tail  velocity.  As 
shown  in  Figure  9,  the  amount  of  the  jet  which  is  characterized  can  significantly  influence  the  measured 
value  of  breakup  time.  The  extremely  high  breakup  time  values  predicted  in  the  tip  region  of  the  jet  are 
due  to  the  mass  buildup,  or  inverse  velocity  gradient,  which  forms  the  tip  particle.  Large  initial  values 
of  jet  breakup  are  measured  experimentally,  although  the  deceleration  of  the  tip  may  be  a  factor. 
Extremely  high  values  of  bieakup  time  arc  also  predicted  for  very  low  jet  tall  velocities  but  these  are 
probably  an  artifact  of  the  Jet  formation  code  since  some  ponion  of  the  liner  (near  the  base)  may  not 
actually  contribute  to  the  Jet. 

Figure  10  shows  the  individual  breakup  time,  calculated  between  the  centers  of  mass  of  adjacent  Jet 
particles,  vs.  the  Jet  velocity.  The  experimental  data  points  show  u  large  amount  of  scatter.  The  scatter 
is  primarily  due  to  variations  in  the  velocity  difference  between  the  adjacent  particles,  as  was  also  shown 
In  Walters  and  Summers  (1992b). 
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Figures.  CumulflUve  breakup  time  as  a  ftinctlon  of  let  temperature  for  the  81 -mm  shaped  charge  with 
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Figure  11.  The  scparuUon  time  of  each  let  particle  aa  a  function  of  let  velocity  for  the  ai-mm 
shaped  charge. 
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Distara^  to  Separation  (mm) 


Figure  12.  The  distance  travelled  before  particulation  vs.  let  velocity  for  the  81-mtn  shaped 
charge. 
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Figures  11  and  12  arc  plots  of  separation  time  vs.  jet  velocity  and  the  distance  a  jet  particle  travels 
before  separation  vs.  jet  velocity,  respectively.  Figure  1 1  indicates  the  actual  particulation  time  tends  to 
increase  with  decreasing  jet  velocity.  It  is  not  evident  that  thw  jet  breakup  proceeds  uniformly  from  the 
tip  to  tail.  The  model  predicts  a  first  break  near  the  front  of  the  jet  but  away  from  the  tip  region  for  the 
warhead  designs  considered  in  this  study.  The  analytical  prediction  agrees  well  with  the  experimental 
data,  however,  the  model  is  clearly  not  sophisticated  enough  to  predict  the  variations  which  occur,  from 
particIc-tO'particlc,  in  individual  breakup  time,  separation  time,  or  separation  distance.  Also,  simple  jet 
formation  codes  arc  not  valid  for  hemispherical  liners,  but  hydrocode  simulations  could  readily  be 
employed  to  determine  the  length  of  various  segments  of  the  jet  at  the  time  they  reach  the  axis  of 
symmetry. 

Figures  13-16  show  a  similar  comparison  between  analytical  and  experimental  data  for  the  140-mm 
conical  liner.  The  agreement  between  calculated  and  experimental  values  is  quite  good  at  the  jet  tail  for 
a  jet  temperature  of  450  K.  However,  the  analytical  predictions  are  not  as  good  near  the  front  of  the  jet 
unless  the  jet  temperature  is  increased  to  650  K.  This  discrepancy  may  Indicate  a  temperature  gradient 
exists  along  the  Jet  (from  the  tip  to  the  tail)  with  the  jet  tip  at  a  higher  temperature.  Table  2  lists  the 
predicted  dp  velocity,  total  Jet  length,  and  cumulative  breakup  time  for  each  of  the  conical  liners  shown 
in  Table  I.  The  total  jet  length  and  cumulative  breakup  time  predictions  are  given  for  the  same 
approximate  velocity  gradient  as  those  given  in  Table  1. 

As  shown  earlier,  the  calculated  breakup  time  results  match  the  experimental  data  quite  well  at  an 
assumed  uniform  temperature  of  450  K.  This  temperature  is  below  the  measured  values  of  von  Holle  and 
Trimble  (Walters  and  Zukas  1989)  for  a  shaped  charge  with  a  copper,  conical  liner.  However,  the 
experimental  temperature  measurements  were  taken  near  the  dp  region  of  the  jet.  There  may  be  a 
temperature  gradient  along  the  length  of  the  jet  (from  dp  to  tail)  as  well  as  a  temperature  gradient  through 
the  jet  thickness.  A  temperature  gradient  from  the  jet  tip  to  the  jet  tall  could  result  in  lower  temperatures 
near  the  rear  ot  the  jet  which  perhaps  undergoes  less  plastic  deformadon  and  less  shock  heating  than  the 
jet  tip  region.  Ideally,  the  individual  breakup  time  could  be  calculated  from  a  jet  formadon  code  with 
temperatures  for  each  element  of  the  jet  (obtained  from  hydrocode  calculations  perhaps)  used  as  input. 

The  breakup  model  in  conjunction  with  the  jet  formation  model  predicts  the  true  strain  at  failure  as 
approximately  2.3.  This  value  is  an  average  value  calculated  for  discrete  jet  segments.  Stress 
concentrations  in  the  necked  regions  of  the  jet  arc  not  included  in  the  model  nor  are  release  waves  or  liner 
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Figure  13.  Cumulative  breakup  time  as  a  function  of  let  velocity  for  the  140>mm  shaped  chfltag 
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Tabic  2.  Jet  Formation  Code  Predictions  for  a  Temperature  of  450  K 


65.3-mm 

Cone 

8 1 -mm 
Cone 

I01.6-mm 

Cone 

140-mm 

Cone 

Vy||)  (km/s) 

9.2 

7.8 

8.9 

8.5 

Vtail* 

4.3 

2.1 

4.2 

2.6 

Lp  (mm) 

508.2 

808.4 

730.3 

1,267.1 

T,,  (MS) 

101.9 

140.4 

154.6 

215.0 

The  tall  velocity  wu  chosen  to  approximate  the  velocity  of  the  last  particle 
characterized  in  the  exporimonti  of  Table  I . 


inhomogenities.  Other  fracture  models  (c.g..  Johnson  [1983])  may  provide  a  more  accurate  strain  at  the 
actual  break  point  between  Jet  particles.  However,  such  a  calculation  would  require  a  hydrocode  solution. 

It  is  also  wonhy  of  mention  that  the  Zciilll-Armstrong  constitutive  equation  has  been  extended  beyond 
its  intended  range  of  applicability  for  the  high  temperatures  and  strain  rates  the  jet  material  undergoes 
during  the  formation  and  particulatlon  process.  Also,  it  has  been  postulated  that  an  increase  in  strain  rate 
is  equivalent  to  a  decrease  In  temperature  (Johnson  and  Mellor  1973).  Thus,  in  the  constitutive  equation 
used  in  this  study,  the  high  strain  rate  effects  may  overshadow  the  temperature  influence. 

5.  DISCUSSION 

As  stated  curlier,  the  average  true  strain  of  the  shaped  charges  with  conical  or  conicul-likc  liners  is 
approximately  2.3.  An  alternative  failure  model  may  simply  be  to  set  the  strain  at  failure  to  2.3.  In  this 
case,  a  priori  knowledge  of  the  Jet  temperature  is  not  required.  Figure  17  is  a  plot  of  the  value  of  true 
stress  predicted  for  a  true  strain  of  2.3  as  a  function  of  jet  temperature  for  both  the  Zerilll>  Armstrong  and 
Johnson-Cook  constitutive  equations.  Other  researchers  have  claimed  failure  occurs  at  a  stress  of  300  MPa 
for  a  copper  Jet  (Walters  and  Summers  1992a).  The  Zcrllli-Armstrong  constitutive  equation  predicts  a 
stress  of  3(X)  MPa  when  the  Jet  temperature  is  approximately  980  K.  The  Johnson-Cook  model  predicts 
a  .stress  of  300  MPa  for  a  Jet  temperature  of  approximately  830  K.  Both  of  these  Jet  temperatures  seem 
high  compared  to  the  measurements  of  von  Hollc  and  Trimble. 
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Alternatively,  the  stress  at  failure  may  be  set  to  300  MPa.  Figure  18  shows  the  values  of  strain  at 
failure  as  a  function  of  Jet  temperature  given  by  the  2^rilll-Armstrong  and  Johnson-Cook  constitutive 
equation:)  for  the  81 -mm  shaped  charge.  Again,  reasonabie  jet  breakup  times  are  achieved  when  the  true 
strain  is  approximately  2.3.  This  faiiure  criteria  gives  reasonable  jet  temperatures,  but  the  strain  at  failure 
is  extremely  sensitive  to  jet  temperature.  A  constant  stress  failure  criterion  reverses  the  trend  shown  in 
the  previous  section  (i.c.,  the  true  strain  at  failure  Increases  with  Increasing  jet  temperature).  In  addition, 
the  trends  predicted  using  the  Johnson-Cook  model  are  in  agreement  with  those  of  the  Zerilii-Armstrong 
constitutive  equation  for  both  the  constant  stress  and  constant  strain  failure  criterions.  For  other  liner 
materials,  a  procedure  similar  to  that  used  in  this  study  can  be  employed  to  obtain  a  new  critical  strain, 
This  strain  value  is  then  constant  for  a  given  liner  material  and  geometry. 

6.  CONCLUSIONS 

A  icchnique  was  derived  to  eoicuiate  the  breakup  time  for  a  shaped  charge  with  a  face-centered-cubic 
liner  materiai.  The  approach  invoived  an  equation  for  plastic  instabiiity;  a  kinematic  relationship  for  the 
jet  breakup  time;  and  u  material-dependent  constitutive  equation,  The  study  was  limited  to  data  which 
the  authors  had  personally  analyzed  since  there  are  several  possible  ways  to  express  a  distribution  of  jet 
breakup  times  and  it  is  often  not  clear  as  to  the  time  origin  Involved  or  the  experimental  data  reduction 
methods  used.  Also,  the  constitutive  equation  used  must  be  expressed  in  terms  of  known  coefficients. 
For  example,  molybdenum  was  not  addressed  in  this  study  since  the  Zerilii-Armstrong  coefficients  are  not 
yet  available  for  this  material. 

The  jet  breakup  time  predictions  were  shown  to  be  very  good  for  several  shaped  charge  designs  with 
copper  line.'}:.  Using  ihc  model  in  conjunction  with  a  jet  formation  code  revealed  good  agreement  between 
calculated  and  experimental  jet  breakup  time  distributions  for  an  assumed  uniform  temperature.  The 
temperature  Is  probably  not  uniform  through  the  jet  but  the  uniform  value  (450  K)  used  in  this  study 
.seems  to  be  low  compared  to  the  measuremenLs  of  von  Holle  and  Trimble. 

The  breakup  time  model  was  further  simplified  by  noting  that  the  average  true  strain  at  failure  of  all 
the  copper  shaped  charge  designs  was  approximately  a  constant,  namely,  2.3.  Also,  a  stress  at  failure  may 
be  selected  (as  300  MPa)  to  yield  another  form  of  the  simpliHed  model. 

Future  studies  will  address  body-centered-cubic  shaped  charge  liner  materials. 


39 


I 


Johnson-Cook 


Zerlltl-Armatrong  - 


Figure  18. 
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